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Abstract 

The electronic and structural properties of pure and defective cerium oxide are investigated 
using a unified LDA+U approach which allows to treat the different valence states of Ce occurring 
for different stoichiometries on a same ground, without any a priori assumptions on the defect 
chemistry. The method correctly predicts the atomistic and electronic structures of the CeC>2 and 
Ce203 bulk phases, as well as the subtle defect chemistry in reduced materials Ce02- X , which 
controls the oxygen storage functionality of cerium-based oxides. The analysis of the energetics 
highlights the limits of the LDA+U method and possible extensions are proposed. 
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Ceria based materials are among the key components of many advanced catalysts used 
to produce hydrogen and to reduce air pollution. The actual catalyst is a complex device 
which couples the high chemical activity of noble metals to the thermo-chemical stability of 
oxide substrates. Besides providing a resistant support for the metal, ceria-based substrates 
take an active role in the catalytic reaction by controlling the oxygen partial pressure at the 
reaction sites, acting effectively as oxygen reservoirs. The accepted qualitative mechanism 
of this phenomenon is based on the following defect chemistry. A reduction of the oxygen 
partial pressure promotes the release from the crystal lattice of atomic oxygen which, after 
diffusing to the surface, recombines and desorbs as O2 gas. The process leaves a charged 
vacancy in the lattice, Vq, which is neutralized by the valence change Ce 4+ — > Ce 3+ of two 
cations (substitutional Ce' Ce defects). The release of one oxygen atom from the substrate 
drives therefore a fundamental change in the electronic structure, leaving two electrons in f- 
like empty states. Electronic correlation, due to the strong localization of these states, have 
the effect of splitting the / band upon occupation, resulting in a fully occupied gap state at 
1.5 eV above the top of the valence band . As a consequence, reduced oxides, Ce02-z, 

all result to be insulating. The electron population of this state is directly correlated to the 
number of Ce 3+ ions, and it is employed as an estimator of the concentration of oxygen 
vacancies 

The localization of electrons on atomic-like / states involves strong correlation effects 
which are not captured by the most standard implementations of density functional theory 
(DFT). As a consequence of this, existing first-principles calculations □ □ □ 3 fail to 
provide a unified picture of the different oxidation states of Ce occurring in ceria. As a 
matter of facts, all of them adopt either one of the following two assumptions, which are 
individually appropriate to different oxidation states of Ce, but conflicting with one another: 
i) in the core-state model (CSM) Ce /-states are treated as part of the core, and therefore 
their contribution to the bonding process is totally neglected j4j ; ii) the valence-band model 
(VBM), instead, treats Ce f electrons explicitly as valence electrons which are therefore 

nriQ 

allowed to contribute to the chemical bond |4|, |y, |2J. The former approach provides good 
structural properties for Ce 2 3 but not for Ce0 2 (see Table|TJ), and it describes the oxide as an 
insulator by construction. Most importantly its predictive power is limited: the distribution 
of Ce 4+ or Ce 3+ ions has to be assumed as an input of the calculation. As a consequence, the 
defect chemistry, which we have seen to be determined by the valence change of Ce atoms, 
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is difficult to access. On the opposite side, the VBM leads to good structural properties for 
Ce02 but not for Ce203 (Table HJ), and describes this latter structure as a metal. Neither of 
the two approaches predicts the existence of the gap state experimentally observed to occur 
in partially reduced ceria CeO^-a, therefore failing to describe the localization-delocalization 
transition of / states, which is at the basis of the oxygen-storage mechanism. 

In the present work these difficulties are overcome using an LDA+U approach which 
allows to describe all the Ce atoms within a single, unified, approximation regardless of 
their oxidation state. This is particularly important in the study of complex, heterogeneous, 
systems — such as e.g. a molecule adsorbed on an oxide-supported metal nano-particle — 
where the oxidation state of each Ce atom is not known in advance. Electron correlations 
within Ce / states are modeled by a Hubbard U contribution added to the traditional 
density-functional self-consistent potential. The method describes explicitly the / states and 
is highly predictive: it correctly describes the electronic structure of pure and defected cerias 
(always predicted to be insulators); it captures the valence change due to non-stoichiometry 
(and the related features in the electronic structure); and it allows the intrinsic defects 
Ce^g to occupy the most energetically favorable position with respect to V' Q without any 
a priori assumptions on their relative geometry. Moreover, it allows to estimate the for- 
mation energetics for pure and defective ceria phases, which turn out to be consistent with 
experiments. 

The results presented in this work were obtained by ab initio DFT-based calculations em- 
ploying the local density approximation for the exchange and correlation potential, and were 
performed with the PWscf package [8]. The crystal valence wave functions were represented 
by a plane-waves basis set limited by the kinetic energy of 30 Ry. Their interactions with nu- 
clei and core electrons were described by non-local ultrasoft psudopotentials [9j constructed 
for the following atomic configurations: O 2s 2 2p A , and Ce 5s 2 5p 6 6s 2 5d 1 4f 1 . The Hubbard 
U contribution to the energy functional was treated after the formulation of Cococcioni and 
de Gironcoli jlOl llll| . The on-site parameter U for Ce was fixed to 3 eV, following the 
analysis of Ref. LL2j where it was estimated with a variational method to lye between 3 and 
3.5 eV. Two concentrations of vacancies were considered in the analysis of defective cerias 
Ce0 2 -:E: x = 0.125 and x = 0.03125. The defects were modeled with the supercell method 
by removing one oxygen atom from the supercells containing, respectively, 24 and 96 atomic 
sites. 
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The LDA+U energy functional reads: 



E LDA+U = E LDA + ^ £ Tr Wi 1 ~ n ')) > (!) 



2 / 

where the n 7 's are M x M matrices (M being the degeneracy of the localized atomic orbital, 
M = 14 in the case of / orbitals), projections of the one-electron density matrix, p, over the 
/ manifold localized at lattice site I: 

(4Lr\p\ < f>L'a') = n mm> S ™'- ( 2 ) 

In the above equation mcr denotes a localized (spin) orbital, with angular and spin quantum 
numbers m and a, and the orthonormality over the spin variables is a consequence of the 
collinearity of the spin structure. The second term in Eq. (fljl, which we call Eu, is positive 
definite for U > because the eigenvalues of the n 1 matrices — i.e. the occupation numbers 
of the / orbitals — lay in the range [0, 1]. Note that, as soon as U is large enough to open a 
gap between occupied and unoccupied / states, the effect of Ejj is to favor the / occupation 
numbers to be close to either or 1. At these extremes, Ejj is strictly and the total energy 
is back to the LDA value. 

Ceria is known to exist in the cubic fluorite-type CeC>2 (Fm3m) and the hexagonal A- 
type Ce2C>3 (P3ml) bulk phases. The latter can be considered as a highly defective CeC>2 
structure where one every four oxygen atoms is missing, with the vacancies being ordered 
along non intersecting (111) directions, and with all the Ce ions nominally in the 3+ valence 



state 5]. The Ce2C>3 structure is anti-ferromagnetic 13], therefore the calculations for this 
phase were carried out in the local-spin-density approximation. 

The calculated structural properties for the bulk polymorphs are summarized in Table |J 
and are in good agreement with experiments. The Hubbard contribution to the energy 
functional increases the lattice parameter by 0.4% in CeC>2 and by 2.5% in Ce2C>3 with 
respect to the LDA values. Note how the LDA+U approach provides a good prediction of 
the structural properties for both polymorphs, therefore combining the achievements of the 
VBM for Ce0 2 with the ones of the CSM for Ce 2 3 . 

The density of electronic states (DOS) for the two polymorphs calculated at the equi- 
librium lattice parameters are shown in Figure ^ Occupied states are indicated by shaded 
areas, and the zero energy is the top of the upper valence band (with prevalent 0-2p char- 
acter). For Ce0 2 (Figure Wfr), the sharp band centered at « 3 eV is formed by fairly 



TABLE I: Calculated structural properties of the bulk Ce02 and Ce203 phases compared with 
previous theoretical values in the valence-band and core-state models (VBM and CSM), and with 
experimental measurements. 





Ce0 2 




Ce 2 3 




a [A] 


B [GPa] 


ao [A] 


c/a 


B [GPa] 


LDA VBM Ref. 5, 


5.39 


214.7 


3.72 


1.55 


208.6 


LDA CSM Ref. 5, 


5.56 


144.9 


3.89 


1.55 


165.8 


LDA+U This work 


5.41 


210.6 


3.81 


1.55 


149.8 


Expt. 


5.41 


204-236 


3.88 


1.55 





localized atomic-like Ce-4/ empty orbitals. The onset of the conduction band is 5.1 eV 
above the highest occupied state, exhibiting the typical LDA underestimate with respect 



to the measured gap of 6.0 eV 



14J (following Ref. [14] we do not consider the /-manifold 



as conduction states). The calculated electronic structure for Ce2C>3 is in good agreement 
with experiments: it predicts an anti-ferromagnetic insulator, with a magnetic moment of 
2.25 /ie/molec (the experimental value is 2.17 /Lie/molec [l|)- The corresponding DOS is 
displayed in Figure [TJd, showing a splitting between occupied and unoccupied states within 
the / manifold. Due to the effects of the Hubbard energy term, the Ce-4/ states do not 
form one single band — as in the CeC>2 structure or in previous VBM calculations — but 
are split in three energy regions: a doubly occupied state at ~ 1.2 eV above the top of 
the valence 0-2p band and corresponding to electrons localized on the Ce 3+ atoms (with 
alternating spins along the z axis). The position and intensity of this peak (relative to the 
0-2p band) are in excellent agreement with experimental photoemission measurements j^J ■ 
In the range of U suggested by Ref. [l^ (3-3.5 eV), the energy of the gap state depends 
linearly on U spanning the values between 1.2 and 0.9 eV (with respect to the top of the 
valence band). The remaining 13 / states are also split in two bands centered at 4.5 and 
7.5 eV, and grouping the 6 and 7 states having their spin parallel and anti-parallel to the 
one of the electron localized on their respective site. 

An oxygen vacancy is formed by removing a neutral oxygen atom from a CeC>2 supercell. 
As a consequence, two electrons have to be accommodated in the / states above the Fermi 
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FIG. 1: Density of electronic states for pure CeC>2 (a) and Ce203 (b) bulk phases, and for defective 
(c) ceria structures Ce02- X - Occupied states are indicated as shaded areas, and the zero energy is 
set to the top of the valence band. 
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energy of Figure Our LDA+U model predicts the localization of two electrons on two 
Ce atoms neighboring the vacancy. This has the effect of formally modifying the valency 
of Ce from +4 to +3. This result is in agreement with the experimental evidence that the 
compensating defects are Ce^g which tend to cluster around the oxygen vacancy, forming a 
defect aggregate. An electronic configuration in which two electrons were initially localized 
on Ce atoms far from the vacancy led to the same self-consistent solution with the Ce 
nearest neighbors to the vacancy being formally Ce 3+ . Upon removal of an oxygen atom, 
the O atoms neighbouring the vacancy relax toward the vacancy by ~ 0.125 A, and the Ce 
atoms outward by ~ 0.08 A. Moreover, the O atoms neighbouring the Ce' Ce relax away from 
these defects (see our previous description on the supercells used to model these systems). 
The DOS for a model of defective ceria, Ce02_ x , is shown in Figure ^ for x = 0.125. 
One sees that the DOS displays in this case features which are intermediate between those 
of the Ce 2 3 and Ce0 2 structures: the gap state at ~ 1.2 eV (due to occupied / states 
on Ce 3+ atoms as in Ce 2 3 ), and the sharp unoccupied / band of the Ce 4+ atoms at 3 
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eV. The position and relative weight of the peak corresponding to the occupied gap state 
are in excellent agreement with observed photoemission spectra . In Figure we 

display the integrated charge density corresponding to the gap state in Ce0 2 _ x . This figure 
demonstrates that, in these dilute regimes, the effects of the formation of a vacancy on the 
electronic structure are limited to the first shell of atoms around the defects (vacancy and 
compensating Ce 3+ ). 

Let us now examine the energetics of reduction of ceria, and focus in particular on the 
two reactions: the Ce0 2 -Ce 2 03 transition, 2Ce0 2 — > Ce 2 3 + l/20 2 (g) (whose enthalpy 



measured at 298 K is —0.58 Ry 15 Q, and the formation of a vacancy (whose experimental 
heat of reduction is 0.34 — 0.37 Ry 16, ^|). The LDA+U approach severely underestimate 
the energetics of these reactions: —0.33 Ry for the Ce 2 3 -Ce0 2 transition, and 0.26 Ry 
for the vacancy formation. This error is not due to the local density approximation itself, 
since gradient corrected calculations (including the Hubbard term) also led to quite similar 
underestimates. 

The source of this error can be traced back to some details of the implementation of the 
LDA+U method which are best evidenced by analyzing the occupation matrices, Eq. (J2j)- 
The integrated charge of the gap state is exactly 2 (one electron for each Ce 3+ atom). 
Inspection of the isosurface charge density displayed in Figures 121 shows that the the gap 
state is not localized exclusively on Ce +3 atoms but has small projections also on the first 
neighbouring O sites. As a consequence, when this state is projected onto Ce-f orbitals, 
Eq. (J2J), the resulting occupation number is clearly smaller than 1 (in fact, ~ 0.98). The 
problem is even more clear in Ce0 2 , where, due to the prevalent ionic bonding, the Ce-4/ 
states should be empty and the occupancies very close to 0. This is not the case, since most 
of the calculated occupancies are small, but not vanishing (in fact,~ 0.01 — 0.02). In the 
particular case of the xyz component, which points along the Ce-0 bond, the occupancy is 
much larger (« 0.2). This large value comes from the partial covalent character of the Ce-0 
bond, has nothing to do with correlation, and should not contribute to the Hubbard energy, 

We conclude that the errors in the vacancy formation energy has to be traced back to 
the spurious occupancies resulting from the particular choice of the projector in ([2|L which 
picks up bonding charge from neighbouring O atoms. Nevertheless, in this relatively simple 
case (where the occupied / state is well separated from the valence band) a simple and 
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effective remedy can be found. If the states used to calculate the occupation matrices, Eq. 
(J2J), instead of being the atomic Ce-f states, were those resulting self- consistently from the 
calculation (i.e. the Wannier functions of prevalent Ce-f character), then the occupancy 
corresponding to the gap state would be one, and that corresponding to other orbitals of 
prevalent Ce-f character zero, by construction. This argument shows that — in the present 
situation where the split state is well separated from all the others — a better characterization 
of the localized orbitals which define Ejj would simply result in a vanishing contribution of 
Ejj to the final, self-consistent, energy. A quick-and-dirt fixup to the vacancy formation 
energy problem seems therefore to be to use Ejj to generate the correct electronic structure, 
but to ignore it altogether in the calculation of the total energy. This amounts to use the 
LDA+U model to generate the molecular orbitals, and revert to the LDA functional for 
the calculation of energies. With this extremely simple scheme, also the energetics of the 
two reactions result in good agreement with experiments: —0.610 Ry for the Ce203-Ce02 



transition (experimental value is -0.58R 
(experimental values are 0.34 — 0.37 Ry 



and 0.361 Ry for the vacancy formation 



We stress that the dependence of the results on the projector used to define the occu- 
pancies is important on the energetics, but have little effect on the electronic structure. 
The atomistic and electronic structures of pure and defective cerias, as well as their ener- 
getics, can therefore be captured by a modified LDA+U method which takes advantage of 
the Hubbard U effect (essential for giving the correct description of the electrons) in the 
self-consistent electronic problem, and whose energy has been "corrected" a posteriori to 
account for the spurious occupancies due to partially covalent character of the chemical 
bond. 

In summary, the LDA+U method provides a simple and accurate way to account for the 
subtle mechanism of /-electron localization-delocalization which determines the mixed va- 
lence character of Ce and which is at the basis of the oxygen storage properties of ceria. This 
method is accurate and predictive, opening the way to the ab-initio study of complex struc- 
tures of technological interest — such as e.g. ceria-supported metal nanoparticles — which 
cannot be easily treated by more empirical approaches where the valence state of indi- 
vidual Ce atoms is pre-assigned. Limitations of this approach show in the calculation of 
reduction/oxidation energies, which can however be fixed in the present case by a simple 
prescription on the energy functional. 
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FIG. 2: Isosurface in the charge density (grey) corresponding to the occupied gap-state in CeOi.969 
(left) and Ce203 (right). Black and light grey circles denote Ce and O atoms. 
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